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Abstract: We discuss the current and prospective limits that can be placed on supersym- 
metric contributions to b — > s transitions in the Minimal Supersymmetric Standard Model 
with general flavour mixing amongst the squarks. We consider three processes: B — * Xg'j, 
B s — ► and B s — B s mixing, and pay particular attention to the large tan/3 regime 

and beyond leading order contributions where the difference between our analysis and pre- 
vious analyses is most pronounced. We find that even present limits on BR(5 S — * 
and AMg s often provide additional constraints on the amount of flavour violation still 
allowed by BR(S — ► -Xyy). Limits on supersymmetric and Higgs mass parameters can 
also be strongly dependent on the flavour violation. In particular, the case \i < can 
still be allowed, even for light sparticle masses. We discuss how future measurements at 
the Tevatron, the LHC and elsewhere can improve present limits but also provide unique 
signatures for the existence of flavour violation. 
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1. Introduction 

The Minimal Super symmetric Standard Model (MSSM) has, for many reasons, become 
one of the most keenly investigated frameworks of particle physics beyond the Standard 
Model (SM). Whilst its ultimate test will be performed at the Large Hadron Collider 
(LHC), it is useful, in the meantime, to place limits on its parameters using existing 
measurements. Examples of such measurements are the W boson mass and sin 2 9 e g, the 
anomalous magnetic moment of the muon and flavour changing neutral current (FCNC) 
processes. 

FCNC processes, in particular, provide a powerful constraint on the MSSM, as the 
SM and supersymmetric (SUSY) contributions each enter at the one-loop level. As the 
SM and SUSY contributions can be comparable to one another, FCNC provide a unique 
means of probing the possible flavour structure of the MSSM. Coupled with the increasingly 
stringent limits that are emerging from the B factories and the Tevatron it is natural to 
investigate how these constraints might affect the available parameter space in the general 
MSSM and how such measurements might complement direct searches for SUSY. In the 
most general formulation of the squark soft terms - namely general flavour mixing (GFM) 
- such measurements can become especially pertinent when constraining the exact form 
the SUSY soft terms can take. 

The usefulness of these constraints, however, is highly dependent on the accuracy of 
the underlying calculation. It is therefore essential to at least to attempt to calculate 
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the SUSY contributions to FCNC processes at a similar level of accuracy as the existing 
SM calculations Q. In general, the additional particle content of the MSSM makes such 
calculations rather prohibitive as one must consider diagrams that, for example, feature 
gluinos as well as gluons, and take care of possible mass hierarchies that exist in the 
particle spectrum, resumming any large logarithms that might occur. Complete next-to- 
leading order (NLO) calculations have therefore only been completed in a few limiting 
cases 0, [|] and tend to focus on the more constrained assumption of minimal flavour 
violation (MFV) [Q, || rather than the GFM framework. With these comments in mind it 
is useful to devise a formalism that builds on existing leading order (LO) calculations to 
include the effects that one might consider to be dominant beyond the leading order (BLO). 
Such effects are typically considered to be proportional to either large logarithms or tan j3 
(the ratio of the vacuum expectation values of the two Higgs doublets that appear in the 
MSSM). Large logarithms are induced by hierarchies that might exist between the coloured 
SUSY particles (the squarks and gluinos) and the characteristic electroweak scale. The 
corrections proportional to tan /3, on the other hand, are induced by threshold corrections 
to the down quark masses (the bottom quark mass in particular) || , and the charged j?], || 
and neutral Higgs vertices pi ITT], IT5] . In the GFM scenario, and without adopting the mass 



insertion approximation (MIA), it has been shown in [13, 14, 15] that the effects induced 
by the inclusion of these tan (3 enhanced threshold corrections can be large, especially 
when compared to similar calculations in the MFV scenario. Broadly speaking, in the 
phenomenologically more viable region fi > 0, the inclusion of such BLO corrections leads 
to a focusing effect that can significantly loosen the bounds on SUSY sources of flavour 
violation [13, [lj|, [TJ|]. For n < the situation is rather more complicated. On the one hand, 
the inclusion of large logarithms can act to decrease certain supersymmetric contributions 
to a given decay (in the case of B — > X s ^ this is particularly true), on the other, the tan/? 
enhanced threshold corrections that acted to suppress the supersymmetric contributions 
for \x > now act to increase them. Broadly speaking, however, it is possible, once one 
includes GFM effects, to allow for /i < 0, even for light sparticle masses (a region of 
parameter space excluded in the MFV scenario), by considering relatively small sources of 



SUSY flavour violation y, |lj. 



The aim of the paper is to present a discussion of the present limits on SUSY flavour 
violation and the future prospects that might arise from measurements at the Tevatron and 
the LHC. We shall pay particular attention to the large tan j3 regime where the differences 



between our analysis and existing LO analyses in the literature [16] are most pronounced. 
The processes we shall use to constrain the various sources of flavour violation in the MSSM 
are: the inclusive decay B — * X s ^; the rare decay B s — * fj> + fi~; and the B s — B s mixing 
system. The SUSY contributions to all three of these processes in the large tan (3 regime 
are known to be sizeable. 



The current world average for the branching ratio B — ► X s j is [17] 



BR(B -> V s7 ) cxp = (3.39±°;30) x 10" 4 . (1.1) 
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When compared with the current SM prediction for the branching ratio [18] 



BR(B -> X s7 ) S M = (3.70 ± 0.30) x 10" 4 , (1.2) 

it is clear that the contributions that arise from SUSY partner exchange in both the MFV 
and GFM scenarios must be small or suppressed in some way. A natural way this can be 
accomplished in the GFM scenario is via the inclusion of BLO effects in the large tan (3 
regime. As discussed in [pH| , 14, 15], when tan/3 is large and /x > 0, sizeable cancellations 



between the chargino and gluino contributions to the decay can occur once one has included 
the various tan/3 enhanced BLO corrections. These cancellations inevitably lead to a 
suppression of the overall super symmetric contribution to the decay and can, in turn, 
significantly loosen the bounds imposed by the decay in the large tan (3 scenario [13, 14, 15]. 



The other two processes we shall consider are currently unobserved. The first is the 
rare decay B s — > / u + fj,~ where only a upper bound exists for the branching ratio. The 
current published 95% confidence limits are JlTJ, |20| 

BR{B S -> /u + /Ocdf < 7 - 5 x 1(r? ' (°) 
BR{B S -» m + aOd0 < 5.0 x 10~ 7 . (1.4) 

Over the past year both CDF and D0 have provided preliminary results for the limit on 



the branching ratio [21, 22] 



BR{B S -> m+zOcdf < 2.0 x 10~ 7 , (1.5) 
BR{B S -» /uV~)d0 < 3.8 x 10~ 7 , (1.6) 



that may be combined to provide the improved limit BR(£> S — > ) exp < 1.5 x 10 |23| 



Bearing in mind that the SM predicts [24] 



BR(B S /u + /i-)sM = (3.46 ± 1.5) X 10~ 9 , (1.7) 

it is clear that the observation of the decay at the Tevatron would provide a clear signal 
for physics beyond the SM. 

The second process we shall consider concerns the mixing between the neutral B® 
mesons. The relevant observable here is the mass difference between the mass eigenstates 
formed between the two mesons AMg s . Currently, only a lower bound exists for the 
quantity [17| 

AM" xp > 14.5 ps-\ (1.8) 



while the SM prediction is |25| 



AMgf = (18.0 ± 3.7) ps -1 . (1.9) 

Comparing the current experimental limit with the SM prediction it is clear that, should 
AA1b s be measured just above the current experimental limit, such values can be typically 
reconciled with the SM prediction ( |1.9[ ). The more phenomenologically interesting region 
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- from the point of view of observing new physics effects - is, therefore, where new physics 
induces contributions to AMb s far in excess of the SM prediction. For example, if AMg s > 
30 ps _1 , such a measurement would imply a deviation from the SM at the level of 2~3cx p(| . 

While only limits currently exist for both of these processes, they already provide a 
useful limit in the large tan/3 regime due to the large enhancement that both processes 
receive. In GFM, in particular, it has been shown that useful constraints can be placed on 
SUSY flavour violation using these processes. 

To summarise what follows, we first describe the bounds imposed by each process in 
turn, discussing the currently allowed regions of parameter space and how future measure- 
ments might affect this picture. In sections ^ and |?] we combine all three constraints to find 
the allowed regions of parameter space for single and multiple sources of flavour violation 
respectively. 



2. Procedure 

Before moving onto our numerical analysis it will be useful to briefly overview the formalism 
that we will use to implement the BLO corrections to the various decays considered in this 
paper, for more details we refer the reader to ]14| , |l5| where the necessary ingredients for 
such calculations are covered in great detail. 

As discussed earlier in this paper, calculations of large effects that appear beyond the 
leading order typically entails the inclusion of the effects enhanced by either large logs 
of the form log M su SY / '^w > where Msusy is a mass scale associated with the coloured 
superpartners and is to be associated with the electroweak scale, or by tan (3. One 
way of including both of these contributions is to consider an effective theory where all the 
SUSY particles have been integrated out. The resulting theory is, essentially, a two Higgs 
doublet model that includes the threshold corrections induced by integrating out the SUSY 
degrees of freedom. To do this it is necessary to assume the following mass hierarchy 

Msusy ~ K, rn~) > fj, w ~ (m H , rn w , m t ) > fi b ~ m b . (2.1) 

That is the squarks and gluino are heavier than the Higgs sector, electroweak gauge bosons 
and top quark, which are in turn heavier than the mass scale of the decays under consid- 
eration. The remaining SUSY particles, such as the neutralinos and charginos, are usually 
assumed to arise in the interval between Msusy and nw 

When working in this effective field theory formalism tan (5 enhanced corrections appear 
as threshold corrections to the masses and couplings that appear in the theory defined at 
Hw As an example, consider the down quark mass matrix. In the effective theory the 
down quark mass matrix m d , which should be identified with the physical mass matrix for 
the down quarks, is related to the bare mass matrix defined in the full theory before 
threshold corrections are taken into account, by the relation 

m d =m^ + 5m d = diag (m d , m s ,m h ) , (2.2) 

where 5m d denotes the threshold corrections calculated in the basis where m d is diagonal 



(the physical super-CKM basis), exact expressions may be found in [15]. The flavour struc- 



ture of 5m d , and therefore m d °\ is highly non-trivial and, even in MFV, the matrices have 
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off-diagonal elements that are enhanced by tan (5, the inclusion of GFM effects only serves 
to enhance these corrections. The non-trivial flavour structure of the bare mass matrix 
acts as an additional source of flavour violation through its appearance in the down squark 
mass matrix and supersymmetric vertices, which are not subject to the threshold correc- 
tions that affect the down quark mass matrix. Similar corrections arise for the couplings 
of the quarks to the gauge and Higgs bosons, the corrections to the Higgs boson couplings 
in particular, can be appreciable even if Msusy is large. The additional contributions to 
flavour violating processes induced by these corrected couplings and masses represent the 
tan/3 enhanced corrections that appear beyond the leading order. It should be noted that 
when one works in the mass eigenstate formalism, where one diagonalises the squark mass 
matrix and flavour violation is mediated by the matrices T^l and r^_R, it is necessary to 
employ an iterative procedure to accurately calculate m?^ and the other vertices in the 



theory, the exact details of this procedure are discussed in |L4|, 15]. 

Now that we have discussed the formalism we shall use let us now discuss the remaining 
aspects of our calculation. After one has calculated the bare mass matrix and corrected 
electroweak vertices the supersymmetric contributions to the process under consideration 
are evaluated, taking into account the effects of the bare mass matrix, and evolved from 
the scale Msusy to the electroweak scale fiw using the relevant six flavour anomalous 
dimension matrix. The remaining electroweak contributions (that is the gauge boson and 
Higgs contributions) are then evaluated using the uncorrected vertices when evaluating the 
NLO corrections and the corrected vertices when evaluating the LO contributions. The 
combined SUSY and electroweak contributions are finally evolved to the scale fif, and used 
to calculate the relevant observable for the process in question. 

Flavour violation in the soft breaking Lagrangian is often parameterised in terms of 
the dimensionless quantities 



■ ■ / j \ ("'-/ , 



SIr).= n I I > ( 2 - 3 ) 



■ "'i/./.).. ('"i/./.) .. ^ 13 \ ( m lLL).A m d,RR, 

m3, o r ) [mh 



, d,RL j . . , , \ \'"'d,RR).. 

3^)^=— ^ J v (8%^^=—= ^ Aj (2.4) 

' fm d,RR) . . ( m d,Ll) . . \l ( m d,RR) . . ( m d,RR) . . 

where i,j = 1, 2, 3 and m 2 d XY (with X, Y = L, R) are related to the soft terms that appear 
in the soft SUSY breaking Lagrangian via unitary transformations that transform the quark 
fields from the interaction basis into the so-called physical super-CKM basis basis where 
the appropriate mass terms are flavour diagonal (for more details see [14, |15|). In the limit 



of MFV all 5xy = 0. As we are primarily concerned with flavour violation between the third 
and second generations we shall, henceforth, use the convenient shorthand 5 ( f jL = (5^) 23 
and so on. 

For the remaining parameters used in our numerical analysis we shall employ a similar 
parameterisation to [15] and treat the soft terms defined in the physical super-CKM basis 
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Figure 1: Plots depicting the effect of varying nig on BK(B — + X s 7) (the left panel) and the 
constraints from BR(_B — > X s 7) on the insertion Sf R (the right panel). In both panels the exper- 
imentally allowed range of BR(_B — ► X s j) is taken at 2a. In the right panel the allowed regions 
corresponding to a BLO calculation are shaded in dark green (dark grey) and bounded by solid 
lines. The equivalent regions corresponding to LO calculation is bounded by dashed lines with 
no additional shading. In both panels the soft sector is parameterised as follows rriq = 1 TeV, 
m A = /i = -A u = 500 GeV and tan = 40. 



as input. For the diagonal elements we set 

[ m 4,LL) u = m \ hi, HUr)*; = m \ 5 Hi { m 4,LIl) u = A d (m d ) u , (2.5) 

while the remaining off-diagonal elements are related to the parameters 5^ Y Yia the rela- 



tions denned in (2.3)-(2.4). The soft terms in the up squark sector are defined analogously. 
As inputs for the Higgs sector we take tua (the mass of the pseudoscalar Higgs), \i and 
tan/3 and use FeynHiggs 2.2 |27j to determine the remaining parameters. For the majority 
of this paper we will only vary one dj^-y & t a time unless stated otherwise. Finally, the 
gaugino soft terms Mi and M2 are related to the gluino mass via the usual unification 
relation. 



3. Limits from B X s ^ 

The first process we shall consider is B — > X s j. Before discussing the impact of the current 
constraint on all four insertions, let us first revisit the focusing effect [jl3| that exists for 
the decay in the large tan /3 regime. Such a situation is illustrated in Fig. |] where we show 
the results of a LO and a BLO analysis beside one another for the insertion <5^ R . As is 
clear from the plots, the difference between the two calculations tends to be rather large. 
In the left panel we show the variation of the branching ratio on m-g for fixed As is 

evident from the figure, for positive and negative 5f R the BLO calculation leads to values 
of the branching ratio that are far more compatible with the current experimental and SM 
results. As discussed in |l3], |l4], this is mainly due to a cancellation that occurs between 
the gluino contribution to the decay (that tends to be reduced by BLO corrections) and 
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Figure 2: Plots depicting the effect of varying fi on the constraint supplied by B — > X s 7. In 
both panels the experimentally allowed range of BR(S — > X s 7) is taken at 2a. In the left panel 
the branching ratio for the decay is plotted against [i. In the right panel, contours corresponding 
to the 2a limits on the branching ratio for BLO and LO calculations are shown. The regions for 
a BLO calculation in agreement (at 2a) with the experimentally allowed range of BR(£? — > X s -f) 
are shaded in dark green (dark grey) and bounded by solid lines. The equivalent regions for a LO 
calculation are bounded by dashed lines with no additional shading when they coincide with the 
allowed region corresponding to a BLO calculation and by light green (light grey) otherwise. In 
both panels the soft sector is parameterised as follows rriq = m g = 1 TeV, vtla = — A u = 500 GeV 
and tan (3 = 40. 



the BLO corrections to the chargino contribution. This cancellation substantially reduces 
the SUSY contribution to the decay and leads to the large reductions in the lower limit on 
Trig evident in the figure. 

The implications the focusing effect has on the limits that can be placed on the insertion 
are illustrated in the right panel, where we show the 2a contours in agreement with the 
experimental result for BR(i? — > X s j) fll.ip 1 in the 5^ R — m-g plane for a LO and a BLO 
analysis. As is evident from the figure the difference between the two calculations can be 
large and the bounds placed on the insertion can be relaxed substantially. 

It is natural to ask how, for fixed nig and m-g, the remaining free parameters affect the 
focusing effect at large tan (5. The largest contribution to the focusing effect usually arises 
from the cancellation that takes place between the gluino and BLO chargino contributions 
to the decay. Regions of parameter space that maximise the chargino contribution to the 
decay with respect to the gluino contribution would therefore be expected to lead to regions 
where the focusing effect is greatest. As the squark masses essentially enter the gluino and 
chargino contributions to the decay, the only difference between the two stems from the 
appearance of the chargino masses rather than the gluino mass in the relevant matching 
conditions. The focusing effect is therefore usually maximised when the two masses that 
enter the diagonal terms of the chargino mass matrix, namely M2 and \x are much smaller 

x To determine the contours we add the SM and experimental errors in quadrature and linearly add an 
additional error of 0.2 x 10 -4 to represent the SUSY aspect of our calculation. 
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than m-g. While the soft mass M2 is typically lighter than the gluino mass in most SUSY 
breaking models, the /i term on the other hand can be much more variable. One should 
note that whilst increasing fi will often decrease the BLO chargino correction, it tends to 
similarly decrease the gluino contribution to the decay through its dependence on the BLO 
correction e s (for the relevant expression see Eq. (4.15) in |15|| ). This compensation tends 
to reduce the overall \i dependence of the focusing effect. 

The \x dependence of the branching ratio is presented in the left panel in Fig. ||. As 
shown in the discussion above, increasing fx tends to have relatively little impact on the 
branching ratio. This is caused by the interplay between the BLO corrections to the gluino 
and chargino contributions mentioned earlier. Similar behaviour is exhibited for RL and 
RR insertions. For LL insertions, on the other hand, the \x dependence is slightly more 
variable as a chargino mediated contribution exists at LO (unlike the other three insertions). 
Turning to the right panel, we show the bounds on the insertion Sf jR and how they vary 
with fi. As is evident from the plot, in contrast to the right panel of Fig. |], the bounds on 
the insertion tend to become more stringent with incerasing jj,. However, we should point 
out that, even in regions of parameter space where /i > nig, the differences between BLO 
and LO calculations can still be large. 

The bounds on each flavour violating parameter imposed by the decay at large tan f3 
are illustrated in Fig. |3[ In the top left panel, depicting the limits on the insertion 5f^ L , 
we can see that the bounds are fairly strict. This is because, at leading order, two large 
contributions to the decay appear. The first is associated with gluino exchange, where the 
dominant contribution stems from the tan/3 enhanced correction that appears at second 
order in the MIA. The second is associated with chargino exchange. As flavour violation 
between left handed squarks in the up and down sectors is related by SU(2) symmetry, 
it is likely that, if is non-zero, the corresponding measure in the up-squark sector - 
5^ L - is non-zero as well. If b\ L is non-zero, it is possible to induce a LO contribution 
via chargino exchange that benefits from appearing at first order in the MIA as well as 
being tan/3 enhanced (for a complete set of analytic expressions see [f!5}j ). The presence 
of this LO chargino contribution tends to negate a large part of the focusing effect as it 
interferes with the cancellation between the BLO part of the chargino contribution, that is 
proportional to 5f L , and the LO gluino contribution. As such, even once BLO corrections 
are taken into account, the impact of the B — ► X s ^y constraint on LL insertions at large 
tan (3 tends to be more stringent than the other three insertions. 

In the top-right panel the bounds on the insertion 5^ R are illustrated. As is evident 
from the plot, the bounds on the insertion are fairly weak. This is primarily due to the 
focusing effect discussed earlier in this section. At LO the constraints on the insertion 
^lr are more strict (see, for example, the plot on the right of Fig. |l]) due to the chiral 
enhancement m,g /nib the leading order contribution receives. However, once one includes 
the dominant effects that appear BLO, cancellations between the LO gluino contribution 
and the BLO corrections to the chargino contribution can lead to a significant relaxation 
of the bounds on the insertion. 

The lower-left and right panels illustrate the bounds imposed on the insertions 5 RL 
and S RR . As is evident from the plots the bounds derived from B — > Xsj on these two 
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Figure 3: Contour plots depicting the impact of BR(£? — > X s j) on the rriq — <5^- y plane. In each 
plot the soft sector is described as follows m g = y/2m,q, A u = —rriq, fx = mq/y2, rriA — 500 GeV 
and tan/3 = 40. The light green (light grey) regions highlight areas that agree within 2a with the 
experimental result, the dark green (dark grey) regions within la. 



insertions are rather slight. This is due to the fact that the contributions due to the two 
insertions cannot interfere with the SM contribution as the dominant contributions due 
to the insertions are to the primed Wilson coefficients. The bounds on each insertion are 
therefore fairly symmetric. For the insertion 5 RL the BLO correction to the charged Higgs 
vertex arising from higgsino exchange leads to slight deviations from this behaviour [15]. 

As an aside let us briefly discuss the decoupled limit m^, m-g — > oo where we are 
effectively left with a two Higgs doublet model that includes the tan (3 enhanced threshold 
corrections discussed in detail in 15]. In this limit the insertions 5f iR and S RL , which scale 
as 1/Msusy j are naturally tiny, and the only meaningful constraint B — * X s ^y can place on 
sources of SUSY flavour violation is on the insertions 5f^ L and o RR . For the insertion 6f L , 
the threshold corrections to the charged Higgs vertex can lead to useful bounds that can 



compliment the constraints supplied by B s — * fi + jjT and B s — B s mixing rather well. For 
the insertion 5rr, on the other hand, the corrections to the charged Higgs vertex only affect 
the primed Wilson coefficients and the constraints that can be imposed on the insertion 
are rather mild to say the least. We shall discuss this issue further at the end of section |6|. 

Finally, let us discuss the scenario fi < 0. Generally in MFV the B — » Y s 7 constraint 
for (j, < tends to be fairly stringent as the chargino and charged Higgs contributions 
interfere constructively with the SM result. In the GFM scenario, on the other hand, 
this situation can be avoided by inducing relatively mild GFM corrections that interfere 
destructively with the other SUSY contributions. For fj, < it is therefore possible to 
reconcile SUSY contributions with experiment for either negative LL insertions or positive 



LR insertions. As the negative \i case was discussed in [13, [T4| (see, for example, Fig. 12 
of |lj]) we shall not cover it in detail here, however, we shall return to it later when we 
supplement the B — > X s ^ constraint with those provided by B s — > fi + ' \i~ and B s — B s 
mixing. 

4. Limits from B s — > /z" 1 " \x~ 

It is well known that SUSY threshold corrections to the neutral Higgs vertex can lead the 



branching ratio for the decay B s — > [J. + [i~ to vary as tan (3 [10, 11, 12 1. Coupled with a 
relatively weak dependence on the SUSY mass scale (the branching ratio varies as \jvn\ 
rather than l/Mg USY ), the constraint supplied by the Tevatron upper bounds on the decay 
can prove to be useful even when the coloured SUSY particles (i.e. the squarks and gluino) 
are heavy. 

Figure f| illustrates the bounds imposed by the decay on each of the insertions for TeV 
scale squark masses and varying pseudoscalar Higgs mass in the large tan (3 regime. All 
four plots illustrate the useful bounds that can be placed on SUSY flavour violation using 
the decay, if ttia is relatively small (compared to the mass of the squarks). 

Before discussing the other aspects of the figures, let us briefly provide a rough recipe 
that allows those investigating flavour violating effects in the squark sector to decide 
whether the B s — > constraint is relevant to their study or not: 

, . . / 40 \ 3 /500 GeV\~ 2 / 40 \ 3 /500 GeV\ -2 
I^lI~0.35 -— ; \5 d LR \ « 0.026 -— ; (4.1) 



, , nn . , 40 \ 3 /500 GeVV 2 . . xd , „ „„ ( 40 \ 3 /500 GeV^ 2 



Let us note that these formulae are highly approximate and serve only as a rough guide 
and should not replace the limits derived from a full numerical analysis. 

The effect of improving the limit on B s — > is also shown in Figure ||. It can be 

seen from all four plots that even the transition from the published D0 result ( |1.4[ ) to the 
preliminary CDF (|F^) result reduces the allowed region in the GFM parameter space by 
a sizeable amount. It is also apparent from the figure that an improvement of the limit to 
roughly 3 x 10~ 8 (a limit achievable at the Tevatron) will provide an excellent constraint on 
GFM in the large tan (5 regime. Finally, let us point out the cancellations that can occur 
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Figure 4: Contour plot depicting the regions of parameter space excluded by increasing the limit 
on BR(B S — > /i + AO f° r tan/3 = 40. The various shaded regions correspond to regions that satisfy 
the following upper limits: 5 x 10~ 7 light green (light grey); 2 x 10~ 7 light-medium green (light- 
medium grey); 3 x 10 -8 dark-medium green (dark-medium grey) and 5 x 1CP 9 dark green (dark 
grey). The soft sector is parameterised as follows m§ = mr g = 1 TeV and /i = — A u = 500 GeV. 



for each insertion. For flavour violation in the LL and LR sectors, direct interference with 
the MFV contributions is possible and values for the branching ratio approaching the SM 
model value are possible for rriA as low as 200 GeV, or less. For flavour violation in the RL 
and RR sectors, on the other hand, the cancellations that can occur are more complicated. 
In general, the insertions act to increase the branching ratio. However, at large tjia regions 
comparable to the SM can become viable. These regions appear due to a cancellation, 
between the Wilson coefficients C' p and Cio, that arises when one calculates the branching 
ratio for the decay. It should be noted that in all four panels the MFV contribution (i.e. all 
= 0) is such that SM-like values of the branching ratio are impossible. This is, of 
course, due to the choice of parameters we make, rather than being a general feature of 
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the MFV contribution at large tan (3. 

The plots in Fig. |] also serve to illustrate how a lower bound on the pseudoscalar 
mass derived from the upper limits on BR(S S — ► // + ^ _ ) is sensitive to the assumed flavour 
structure, and how any upper limit that appears (should B s — ► fi + ^~ be observed at the 
Tevatron for example) might be altered beyond MFV. As an example let us consider the 
top-left and bottom-right panels in Fig. The top-left panel illustrates that one can 
typically avoid the lower bound on itia altogether even for a relatively minor amount of 
flavour violation in the LL sector. While this assumption smacks of fine-tuning, one should 
remember that LL insertions are by far the easiest to generate through RG running in the 
MSSM (see, for example |28|]), and such notions, therefore, might not be as distasteful as one 
might initially think. Beyond these cancellations that appear for small Sir it is apparent 
from both panels that the presence of flavour mixings in the squark sector generally act to 
increase the lower bound on itia- 



It was discussed in [29] that if one assumed MFV and imposed an upper bound on 
tan/3, it would be possible to place some sort of upper bound on the pseudoscalar Higgs 
mass niA- If one allows for GFM in the squark sector it is apparent from all four plots 
that this result will no longer hold (one can travel along a given band in the plots towards 
arbitrarily large m^) and one would probably need to resort to a more general analysis 
that takes into account a variety of other processes before any concrete bound could be 
derived. 

5. Limits from B s — B s Mixing 

Large contributions to B s —* /i + /i~ have been known for some time now to induce similar 
corrections in the B s — B s mixing system. These corrections arise through Feynman dia- 
grams that feature two of the flavour changing neutral Higgs penguin vertices that feature 



in the tan/3 enhanced contributions to B s — > \v \i~ (see [15[ for more details). In the 
MFV scenario such corrections typically lead to reductions to AM Rs taking it closer to 
the current experimental limit [Bol. Since LL insertions affect the corrected neutral Higgs 



vertex in a similar manner to the MFV corrections [15] the overall effect is similar to that 
found in MFV. The only exception is that the effects are generally more exaggerated due to 
the enhancement that the GFM correction receives due to the presence of the strong cou- 
pling constant in the relevant matching conditions. On the other hand, the contributions 
due to the insertion 5f iR tend to be small and generally are bounded by the B s — ► / u + fi~ 
constraint [Ol. 



Turning to flavour violation in the RL and RR sectors, as we pointed out in [15], the 
contribution where one neutral Higgs penguin is mediated by chargino exchange and the 
other by gluino exchange (thereby eliminating the suppression by m s that blights the MFV, 
LL and LR contributions), can yield large effects on AM Ra - Such effects are illustrated 
in Fig. H where we show contours of AMg s for varying S RL , 5 RR and tua for large tan/3. 
As is evident from both panels there is far more variation in A M Ra compared to the MFV 
limit = 0- Negative values of 5 RR and positive values 6 RL , in particular, seem to 
be disfavoured by the current lower bound on AM Rs - It should be noted that the effects 
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Figure 5: Plots illustrating the dependence on tua and Sx Y 01 AMb s once the published D0 limit 
BR(_B S — > < 5.0 x 10~ 7 has been taken into account (the region shaded in orange/medium 

grey and bounded by a thick line). The shading is as follows: white regions are excluded by the 
lower limit AMb s > 14.5 ps _1 ; light green (light grey) regions indicate values of AMg, in the 
interval 14.5-22 ps _1 ; medium green (medium grey) 22-30 ps -1 ; dark green regions (dark grey) 
indicate points in parameter space where AMb, > 30 ps -1 . The soft sector is parameterised in the 
same manner as Fig. 0. 



induced by varying arise entirely from BLO corrections and would be absent in a 
purely LO calculation Jlq]. From the plots it should also be noted that values for AMb s 



far in excess of the SM prediction (1.9) are also allowed by the current upper bound on 
BR(i? s — > [j, + fi~Y This provides a unique signature for the existence of flavour violation 
in the RR and RL sectors if tan f3 is large. For example, if the Tevatron were to measure 
BR(i? s — > fi + at a level far in excess of the SM prediction, and AMb s was subsequently 
measured at a value substantially larger than the SM prediction (rather than lower as 
predicted by MFV |3(|), such observations could only point towards the existence of non- 
minimal flavour structure in either the RL or RR sectors. 



It has been known for some time [11, 15] that if flavour mixings appear in both the LL 



and RR sectors then the contributions to AMg s can be especially large. Such a situation 
is depicted in the left panel in Fig. ||. Here we can see regions where AMb s exceeds values 
of 250 ps -1 (a value unobservable at LHCb) and still satisfies the current published D0 
bound BR(S S — > fi + fi~) < 5.0 x 10 -7 . Similar contributions arise if flavour mixings appear 
in both the LR and RL sectors as illustrated in the right panel in the figure. Here the 
constraints and the general behaviour of the panel arise only when one includes effects 
that appear beyond the leading order. (One should note, however, that we do not include 
the constraint supplied by B — > X s -f in both panels.) The large contributions that appear 
in both panels only arise if flavour violation in the LL or LR sectors occur together with 
flavour violation in the RL or RR sectors. 
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Figure 6: Plots illustrating the variation of AMb s once the published D0 limit BR(B S — > < 
5.0 x 10 -7 has been taken into account (the region shaded in orange/medium grey and bounded by a 
thick line). The panels are shaded as follows: white regions are excluded by the lower limit AMb s > 
14.5 ps _1 ; light green (light grey) regions indicate values of AMg a in the interval 14.5 ps _1 -22 ps _1 ; 
medium green (medium grey) 22 ps _1 -250 ps _1 ; and dark green (dark grey) regions indicate points 
in paramter space where AMb s > 250 ps -1 . The soft sector is parameterised in the same manner 
as Fig. [|. 



6. Limits on Single Sources of Flavour Violation 

With the limits that can be provided by each process in mind, let us now combine all three 
constraints discussed above. Such a situation is shown in Fig. |7] where we vary m-q along 
with a single insertion. Unlike the figures shown in the previous three sections we use white 
to depict allowed regions in parameter space (rather than green/grey) and choose instead 
to shade the regions excluded by each successive constraint we impose. 

The top-left panel in the figure illustrates the allowed parameter space for the insertion 
It is evident from the figure that the constraints supplied by B s — » fJ, + fi~ and B s — B s 
mixing play a useful role in constraining large values of 5f L . In particular, they tend to rule 
out the extreme regions of parameter space, where the SUSY contributions to B — > X s ^f 
have effectively flipped the sign of the underlying amplitude relative to the SM result. It 
is also apparent that, when one inspects the contours that outline the regions allowed by 
these constraints, both processes display a weak dependence on m,. This is due to the 
decoupling effect of the SUSY corrections to the neutral Higgs vertex [10, 11, 12]. These 



constraints will therefore remain even when rrig approaches values for in excess of 1 TeV 
(provided that tua remains constant). 

The constraints on the insertion 5f R are illustrated in the top-right panel. In this 
plot we see the constraint supplied by B s — > fJ, + [J,~ once again tends to disfavour larger 
values of Sf^ R . It is also apparent from the plot that the constraint seems to become more 
severe with increasing m<r. This behaviour can be readily understood by recalling that the 
insertion Sf jR scales as l/MgusY- If one were to take into account this behaviour somehow, 
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Figure 7: Contour plots showing the current limits on the GFM parameters S^y f° r varying rriq 
and fixed rriA- The soft sector is parameterised in a similar manner to Fig. ||. Regions excluded by 
B — » X s 7 are shaded in yellow (light grey), the subsequent regions that are excluded by B s — > 
are shaded in orange (medium grey), finally the remaining regions that are excluded by the limit 
on AMb s are shaded in red (dark grey). 



by plotting the constraints on the off-diagonal elements of the actual trilinear soft terms, 
for example, the correct decoupling behaviour would be restored. 

The bottom-left plot in Fig. ^ illustrates the constraints on the insertion Sj^j. Here 
we see that the processes B s — ► /U + /i~ and B s — B s remove additional ranges of <5^ L 
for rriq as small as 550 GeV and determine the limits on the insertion completely above 
rriq ~ 700 GeV. Once again the limits imposed by these constraints increase with rriq, 
however, in a similar manner to the insertion 5f iR discussed above, it is apparent that, once 
one takes into account the dependence of the insertion on the underlying SUSY mass scale, 
the correct decoupling will be restored. 

Finally, consider the bottom-right panel in the figure that illustrates the bounds on the 
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sector is parameterised in a similar manner to Fig. | with m A = 500 GeV. The excluded/allowed 
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insertion 5^^. Here we see that, apart from rriq ~ 500 GeV - where MFV is disfavoured by 
the B — ► Xg'j constraint, the limits on the insertion are determined completely by the limits 
on BR (B s — ► fi + 1^ ) an d AMb s - In addition, in a similar manner to the insertion 5f L the 
dependence of these two constraints on the underlying SUSY mass scale is extremely weak. 
As such, the constraints supplied by these two processes will continue to be useful, even if 
the squarks and gluino effectively decouple. 

The SUSY contributions to all three processes under investigation in this paper are 
well known to be highly dependent on tan/3. The corrections to BR (i? s — > fi + fi~^ and 
AMb 3 for instance scale as tan 6 and tan 4 /? while the BLO focusing effect, discussed 
in []l3|, 0, |l5| and at the beginning of section is also strongly dependent on tan/3. It is, 
therefore, natural to ask how the limits on the various sources of flavour violation depend 



on tan (3. Such a situation is shown in Fig. || where we illustrate the allowed regions in the 
parameter space formed by varying tan (3 and one insertion at a time. 

In the top-left panel the constraints on the insertion 5^ L are shown. Here we see that, 
in contrast to the other three panels in the figure, the constraints on the insertion increase 
dramatically with tan [3. For example, at low tan (3, the insertion is relatively unconstrained 
and values of up to 5f L ~ 0.8 are easily possible. This is because the dominant LO chargino 
and gluino contributions to the decay that arise from LL insertions are both proportional to 
tan (3. As tan (3 increases the constraint imposed by B — ► X s j therefore becomes much more 
stringent and separates into two distinct branches. The upper branch here is associated 
with regions where the amplitude associated with the decay has effectively changed sign. 
Once we include the constraints supplied by B s — > and B s — B s mixing, however, 

these extreme regions are typically ruled out for tan [3 > 30 (and also, to some extent, for 
tan/3 < 15). 

The limits on the insertion 5^ R are shown in the top-right panel. Here we see the 
focusing effect rather clearly: the regions of parameter space excluded by the B — > X s ^ 
constraint decrease dramatically as tan (3 increases, in contrast to a LO analysis where, 
essentially, the bounds on are relatively independent of tan [3. As tan (3 increases, 
however, the constraint supplied by the decay B s —* also gains prominence and 

becomes important in determining the limits on the insertion above tan/3 ~ 30. Indeed, 
above tan/3 ~ 55, the constraints on 5f^ R are determined entirely by the published D0 
bound BR(5 S -» fJ, + fJ,~) < 5.0 x 10~ 7 . 

The lower-left panel in the figure illustrates the bounds on the flavour violation in the 
RL sector. Once again, the bounds on the insertion supplied by the B —* X s ^ constraint 
decrease steadily as tan/3 increases. This is due to the focusing effect discussed at the 
beginning of section ||. It is also apparent from the figure that the constraints supplied 
by B s — > fJ. + /J~ and B s — B s mixing play large roles when determining the limits on 5 RL 
at large tan/3. The limits imposed by the current lower bound on AMb s , for instance, 
start eliminating large regions of parameter space for positive 5 RL for tan (3 as low as 30. 
Increasing tan [3 only serves to increase the effectiveness of the constraints supplied by the 
decay B s — > fi + and B s — B s mixing, and, for example, above tan/3 ~ 45 the constraint 
provided by B — > X s ^ is entirely supplanted by those attributable to B s — > and 
B s — B s mixing. 

Finally, in the bottom-right panel we show the limits on the insertion 5 RR . From the 
plot it is evident that the bounds imposed by B — > X s ^f are fairly weak. This is principally 
due to two reasons: the first is that the insertion S RR only contributes (in a significant 
manner) to the primed Wilson coefficients and, therefore, cannot interfere directly with 
the SM contribution; the second is that the BLO focusing effect significantly reduces the 
tan/3 enhanced LO contribution that arises from gluino exchange. It is therefore apparent 
from the figure that the constraints supplied by both B s — > and B s — B s can become 

important even for tan/3 as low as 20 and for tan/3 > 30 the constraints on the insertion 
are determined entirely by these two processes. 

To illustrate the usefulness of combining the constraints supplied by B — > X s ^, B s — > 
and B s — B s mixing we provide a list of the bounds on each insertion in table [l]. 
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s d LL (w- 2 ) 


tan/3 = 10 


tan /3 = 40 


tan/3 = 50 


A 


(-2.2,11.0), (47,56) 


(-4,11) 


[-1.2,-0.2) 


B 


(-10,40) 


(-7,10) 


(-7.4,8.0) 


C 


(-5,8.6), (45.6,54.8) 


(-4.4,0.2), (14.4,18.8) 


(-4.6,-0.7), (12.0,15.8) 


D 


(-17.4,37.0) 


(-9.2,9.6) 


(-8.8,7.4) 




tan/3 = 10 


tan = 40 


tan/3 = 50 


A 


(-1.2,7.6), (22,31) 


(-12,2) 


[-9.2,-1.6) 


B 


(-3,15), (46,64) 


(-12,18) 


[-14.8,9.2] 


C 


(-3.0,5.4), (21.0,29.6) 


(-13.4,0.2), (26.4,33.6] 


(-19.2,-3.6) 


D 


(-5.00,12.8), (45.0,63.0) 


(-14.8,19.6), (55.2,78.8] 


(-20.2,24.4) 


4l(1(T 3 ) 


tan/3 = 10 


tan = 40 


tan/3 = 50 


A 


(-7.0,7.0) 


[-12,0.5}, {6,10] 


Excl. 


B 


(-16,16) 


[-30,2}, {14,26] 


[-16.0,0.6}, {6.6,14.0] 


C 


(-9.8,9.8) 


(-32.8,3.4}, {26.2,32.0) 


[-28.0,-15.0), (13.8,24.8] 


D 


(-19.4,19.4) 


(-52.8,7.6} 


[-64.4,3.4}, {27.6,57.2] 




tan/3 = 10 


tan 13 = 40 


tan/3 = 50 


A 


(-24.6,26.8) 


[-5,-2}, {-0.4,6] 


Excl. 


B 


(-80,84) 


[-23.2,-11.4}, {-1.6,26.4] 


[-13.0,-6.0}, {-0.6,14.8] 


C 


(-36.4, 38.6) 


[-22.6,-11.6}, {-1.6,25.8] 


[-12.6,-6.1}, (9.4,14.6] 


D 


{91.2,92.0} 


(-71.4,-39.2}, {-6.4, 74.2) 


[-50.8,-21.8}, {-3.1,56.8] 



Table 1: Table of limits on all four insertions for varying tan/3. The soft sector is parameterised 
as follows: in all four scenarios rriq = rrig = fj, and A u = —rriq; in scenarios A and C rriq = 500 GeV; 
in scenarios B and D rriq = 1 TeV. The pseudoscalar mass itia is given by itia = rriq/ 2 in scenarios 
A and B, and by rriA — rriq for scenarios C and D. The source of each bound is indicated by varying 
styles of bracket, conventional round brackets indicate B — ► X s j, square brackets B s — > fi + fi~ and 
braces indicate B s — B s mixing. 

The first column in the table illustrates the bounds on the insertions for tan (3 = 10. It is 
evident here that, for this moderate value of tan/3, the constraints on all four insertions 
are determined entirely by the decay B — ► X s ~f (except for one case for 5 RR ). This is, of 
course, to be expected due to the large dependence the neutral Higgs penguin contributions 
to B s — ► and AMb s have on tan/3. Moving on to the next two columns which 

display the limits on the insertions for much larger values of tan [3 it can be seen that 
both B s — > /i + /U™ and B s — B s mixing play a much more useful role. For example, for 
LL and LR insertions the constraints either remove the second set of allowed values that 
were present at tan/3 = 10 or reduce them somewhat, even for relatively large txia- For 
RL and RR insertions the constraints are determined almost entirely (with one or two 
exceptions) by the current limits on BR {B s — > /U + /i~) and AMb s . The allowed values of 
these insertions also take a different form compared with those obtained for small tan/3. 
For low tan (3 the allowed regions are usually given by (roughly) symmetric regions centred 
on the MFV limit (i.e. Sj^y = 0)- However, at large tan/3 the allowed regions appear to 
split into two separate branches. These branches that appear at negative ^^.r positive 
are, essentially, regions of parameter space where the SUSY contributions have flipped 
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Figure 9: Contour plots showing the current limits on the GFM parameters 8f L and Sf jR for 
varying niq and fj, < 0. The soft sector is parameterised in a similar manner to Fig. | while the 
shading of the figure is that same as Fig. [?]. 



the sign of the underlying amplitude, relevant to B s — B s mixing, relative to the SM result. 

Let us now consider the case ^ < 0. As discussed earlier in this paper, for /i < 0, the 
chargino and charged Higgs contributions to B — > X s ~f interfere constructively with the 
SM result and one is often forced to adopt a very heavy SUSY mass spectrum {Msusy 3> 
1 TeV) to ensure that the contributions to BR(£? — > X s j) are not too large. In GFM 
it is possible to improve the situation by considering contributions from either LL or LR 
insertions that interfere destructively with the chargino and charged Higgs corrections to 
the decay. This problem is made somewhat easier when one considers BLO effects as, 
instead of reducing the GFM contribution, as they do for positive fi, they act to increase 
it (i.e. "anti-focusing"). As such it is possible that relatively small deviations from MFV 
will satisfy the B — > X s ^f constraint. As changing the sign of \i also changes the sign of 
the term e s that appears in the various BLO corrections presented in [15|, it is apparent 
that the contributions to both B s — > [i^ [i~ and B a — B s mixing will increase as well. It is 
therefore to be expected that the constraints imposed by these two processes will play an 
even greater role compared to the positive /x case. Such a situation is shown in Fig. ^ where 
we present the allowed parameter space in the 5^-rriq and S^-rriq planes for negative \x. 
Each plot displays the familiar structure of two branches that are compatible with the 
B —* X s ^/ constraint encountered earlier in Fig. ||. In a similar manner to the plots found 
in Fig. |7|the constraints supplied by the limits on BR(I? S — > fi + [J.~) and AMg s tend to rule 
out the more extreme regions of parameter space. For example in each figure the constraint 
supplied B s — > excludes the extreme branch of the region permitted by B — > X s ^y. 

In the left panel of Fig. || the constraint supplied by B s — B s mixing also plays a large role 
and rules out over half of the allowed parameter space in the upper branch in the figure. It 
is obvious from both figures however that regions consistent with all three constraints are 
still possible and fi < can still be a viable possibility at least when one considers minor 
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Figure 10: Contour plots depicting the limits on the insertion 8^ L in the decoupled limit rriq ^> mj. 
The soft sector is parameterised as follows rriq = = rrig/^/2 = —A u = lOTeV and tan/3 = 40. 
In the left panel n > 0, while in the right panel fj, < 0. 



departures from MFV. 

Finally, let us consider the decoupled limit mq,mg 3> tha- Such a situation is illus- 
trated in Fig. 1C where we depict the constraints that can be placed on the insertion 5^ L . 
As is evident from the two panels all three processes continue to play a role in constrain- 
ing the LL insertions. It should be noted that the constraints arising from B — > X s ~f in 
particular, are a purely BLO effect. The constraints arise from the threshold corrections 
to the charged Higgs vertex and can act to either increase or decrease the contribution to 
B — > X s "f arising from charged Higgs exchange relative to a purely LO calculation. This 
is in contrast to an MFV calculation where the BLO corrections act only to decrease the 
charged Higgs contribution for \i > and increase it for \i < 0. It is also evident from 
both panels that the constraints provided by B s — > (J, + [J,~ and AMb s , which arise from 
the threshold corrections to the neutral Higgs vertex, continue to play an important role 
eliminating roughly symmetric regions of parameter space in each panel. Finally, let us 
note that in each panel the lower bound on tha in the MFV limit is determined solely by 
B — > Xsj, however, once one proceeds beyond this approximation the lower limit on tua 
can be dictated by either the B s — * (JT or B s — B s constraints (especially when \i > 0, 
as in the left panel). 

It is natural to ask how the other insertions are constrained in this extreme region 
of parameter space. As the threshold corrections to the charged Higgs vertex, arising 
from RR insertions, only effect the primed Wilson coefficients [|15| it is apparent that the 
constraint supplied by B — > X s -y will have a relatively minor dependence on 5j^ R . As such 
the constraints on the insertion are usually determined by B s — > /U + /i~ and B s — B s mixing, 
and eliminate regions of parameter space similar to those found in Figs. |[-||. Turning to 
flavour violation in the LR and RL sectors, as the original definitions for these scale as 
1/MsuSYi it is generally more useful to constrain the actual elements of the trilinear soft 
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terms themselves, or at least a dimensionless quantity that does not feature a Higgs VEV 
or quark mass. Provided one takes this into account it is apparent from the decoupling 
behaviour exhibited by the B s — ► and B s — B s constraints in Fig. ^ that both 

processes will continue to be useful when constraining flavour violation in the LR and 
RL sectors, even for rriq = lOTeV. Turning to the B —* X s ^ constraint, neither insertion 
induces large threshold corrections to the charged Higgs vertex and as such the dependence 
of the B — > Xgj constraint on these two insertions would be expected to be relatively small 
when compared to the large dependence on 6f^ L exhibited in Fig. [l0|. 



7. Limits on Multiple Sources of Flavour Violation 

Up until now, for the sake of simplicity, we have been mainly concerned with one source of 
flavour violation being varied at a time. However, it is more natural to expect more than 
one insertion to be simultaneously non-zero. This may lead to a more varied picture. For 
example, as discussed at the end of section ||, large contributions to AM b s are possible if 
two insertions are both non-zero. Coupled with the possibility that interference between 
the various contributions can also play a role for B — * X s ^y and B s — > it is interesting 

to consider the limits attainable if one decides to vary two sources of flavour violation in 



the squark sector at once. Such a situation is shown in Fig. 

First let us consider the regions ruled out by the B — ► Xg'j constraint (the yellow 
regions). The regions permitted by the constraint are typically circular in nature with the 
exception of the 6 RIy — S RR plane (the bottom-right panel). The reason for the deviation 
is due to the fact that neither insertion contributes significantly to the unprimed Wilson 
coefficients CV and C% and the overall effect of the insertions is to increase the branching 
ratio for the decay. For the remaining combinations at least one of the insertions contributes 
to the unprimed sector and can therefore lead to a decrease as well as increase in BR(£? — > 
Xsj). Let us briefly also mention why some of the circular regions are filled (mainly those 
featuring the insertion S^r) w hilst others are unfilled (those that feature the insertion 5^). 
The reason here is due to the focusing effect discussed at the beginning of section ||. In 
a similar manner to the top-right panel in Fig. |||, the focusing effect raises the minimum 
value for BR(i? — > X s j) in the circular regions towards the SM value and above the lower 
limit we impose (a more stringent lower limit would lead to "doughnut" shaped regions 
like the remaining panels). As the contributions due to LL insertions are less affected by 
BLO corrections (compared with the other three insertions) the allowed regions in the plots 
featuring this insertion tend to retain a similar form to their LO counterparts. 

All six panels in the figure easily illustrate the usefulness of the current bound in 
B s — » /i + /U _ in constraining GFM at large tan/3 (the dashed lines). It is also apparent that 
the regions excluded by the constraint are rather different in the Sf iL Sf jR and the S RL -5 RR 
planes (the top-left and bottom-right panels respectively) compared with the remaining 
four combinations. This is because, in both cases, the contributions to the neutral Higgs 
penguin interfere directly with one another, leading to the rectangular shaped regions in 
the two panels. It should also be noted that these regions are orientated at roughly ninety 
degrees to the B — * X s ^ constraint in both panels. As such, the combination of the two 
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Figure 11: Shading the same as Fig. 0, the soft sector is parameterised in a similar manner to Fig. 
with m A = 500 GeV. 



constraints tends to be fairly complementary when constraining these two combinations. 
For the remaining four contributions, one insertion affects the right -handed Higgs vertex 
while the other affects the left-handed vertex. These contributions translate to contribu- 
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tions to the unprimed and primed Wilson coefficients respectively. Direct interference is 
therefore ruled out, resulting in the circular regions found in the remaining four plots. 

Let us now consider the constraint imposed by the B s — B s mixing system. Here we 
see that, in all six panels, the constraint imposed by this process is comparatively mild 
(with the exception of the plane). This is mainly due to the way we impose the 

constraint. Since we simply require that AM Bs > 14.5 ps _1 the regions excluded in the 
panels involving the combination of an LL or an LR insertion, and a RL or a RR insertion, 
tend to be rather slight due to the large effects that are possible for these combinations 
(see the end of section ||). Therefore, if one chose to impose some sort of upper bound on 
AMb s as well, the constraints on SUSY flavour violation imposed by the B s — B s mixing 
system would be far more stringent. 



Finally, let us briefly make one more remark concerning how the plots in Fig. 11 



might change should one choose to vary other sources of flavour violation in addition 
to those already being varied in each plot. Generally, it would be expected that quite 
large deviations from the contours shown in Fig. |ll| occur once one proceeds beyond this 
limit, however, there is one exception. Certain regions in the top-left panel of Fig. |ll] are 
relatively independent of the effects of varying <5^ L and <5^r. This is due to the fact that 
varying <5^ L and 5^ R invariably only affects the primed Wilson coefficients. The result of 
varying these insertions is that they, therefore, tend to increase the values of BR(i? — * Xgj) 
and BK(B S — ► fj. + n~) compared to a calculation where they are set to zero. Varying these 
insertions might therefore fill in the excluded region at the centre of the "doughnut" shown 
in the plot. The variation, however, will not affect (in a large manner) the regions excluded 
by the outer ring that delineates the B — » X s j constraint nor the regions excluded by the 
B s — ► constraint. The regions excluded by the lower bound on AMb s , however, 

will tend to vary much more as varying RL and RR insertions together with LL and LR 
insertions inevitably leads to large deviations in AMg s . 

The correlation between certain observables in the large tan (3 regime is well known. 
In |3(J, for instance, a strong correlation between BR(1? S — > /i + /x~) and AMb s was pointed 
out that exists due to the dependence both processes have on the neutral Higgs penguin. 
It was shown in [15] that similar correlations exist if one varies only one insertion at a time 
in the GFM scenario. Fig. 12 illustrates the result of varying all four insertions at the same 
time, in a particular region of parameter space. As is evident from the figure, the correlation 
between BR(.B S — > and AMg s that exists in MFV and certain limits of GFM is less 

pronounced once one varies all four sources of flavour violation. This is hardly surprising 
once one notices the vastly different behaviour exhibited by the contributions arising from 
each insertion to AM Bs [fl"5| l, as well as the fact that multiple sources of flavour violation 
lead to exceptionally large contributions to AMb s (see the end of section ||). The scatter 
plots in Fig. [l^ also demonstrate the effect that increasing the bounds on BR(i? s — > /i + /i~) 
or AMb s might have on the available parameter space. In the left panel (corresponding 
to a relatively light SUSY spectrum and Higgs sector) it is evident that the majority of 
the points consistent with B — * X s ^f have already been ruled out by the Tevatron limits 
on BR(i? s — > Turning to the right panel, where the mass spectrum is twice as 

heavy, it is evident that there is far more freedom, this is of course unsurprising as the 
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Figure 12: Scatter plots illustrating the relationship between BR(i3 s — > and AMfl a when 

all four sources of flavour violation are being varied. In each plot and are varied over 
the range [—0.8,0.8], while the insertions (5£ fl and <5^ L are varied over the range [—0.08,0.08]. 
The soft sector is parameterised using the relations rriq = —A u = mg/\/2 = \/2 (i — 2 m a with 
TYiq — 500 GeV in the left panel and trig = 1 TeV in right panel, tan/3 = 40 in both plots. Blue 
(grey) squares denote points that are consistent with the experimental result for BR(£? — * X s j), 
while black dots are not. 



contributions to all three processes decrease as the mass spectrum increases. It is evident 
from the plot, however, that reducing the limit on BR(S S — ► /x + /i~) to roughly 3 x 1CP 8 
or observing AMb s will affect the available parameter space dramatically. 



8. Other Constraints 



Before proceeding with our conclusions let us briefly discuss how some additional con- 
straints might affect the available parameter space in the GFM scenario. 

The first process we shall consider is the decay B — ► X s l + l~. The effect this process 
might have on the SUSY parameter space was discussed in [[H]] where it was shown that 
the combined constraints provided by B — > X s l + l~ and B — > X s ^f indicate that the sign 
of the amplitude for the decay B — > X s j was that of the SM, unless large contributions 
to the Wilson coefficients Cg S and C^q could be induced. Recalling the discussions in 
section |3] and section |(| it is possible for both LL and LR insertions to induce contributions 
to B — ► X s ^y that flip the sign of the underlying amplitude (see, for example, the upper 
branches of the top-left and top-right panels in Fig. ||). As we have seen, the constraints 
imposed by B s — > /^ + /U~ and B s — B s mixing typically exclude the majority of these regions 
at large tan (3. However it is still possible, in some cases, that such regions can remain even 
after one has taken into account both of these constraints (see, for example, the top-left 
panel in Fig. 11). In addition, it is interesting to consider how one might remove these 
regions at low tan/3 where the B s — > and B s — B s mixing constraints play only a 
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minor role. While a complete analysis of the B — > X s l + l~ constraint in GFM is beyond 
the scope of this analysis, let us briefly make a few comments on the possible effects of the 
constraint. 

For LR insertions it is difficult to induce large contributions to either Cg S or Cfg . This 
is because the operators associated with these Wilson coefficients feature the Dirac bilinear 
SLjubL- The contributions to these operators arising from LR insertions, therefore, appear 
at second order in the mass insertion approximation and suffer from a suppression by either 
the strange or bottom quark mass. Compared with the large contributions induced by LR 
insertions to the Wilson coefficient C| ff that appear at first order in the MIA, in addition 
to being chirally enhanced by the gluino mass, it is likely that the constraint imposed by 
B — > X s l + l~ will be useful in constraining the large values of LR insertions that have 
effectively flipped the sign of the B — > X s ^ amplitude, at low tan j3 in particular. 

Turning to LL insertions, it is possible that contributions to the operators associated 
with the Wilson coefficients C| ff and C^q might be somewhat larger as contributions pro- 
portional to Sj^l can appear at first order in the MIA via photon penguins [32]. These 
contributions are not enhanced by tan /3, however, and it seems likely that for large tan f3 
the constraint will play a similar role to B s — ► and B s — B s mixing, due to the tan (3 

enhancement of the corrections to C| . The impact the constraint might have on con- 
straining RL and RR is less clear as they affect the primed coefficients and the possibility 
of flipping the sign of the B — ► X s ~f amplitude is removed. 

In summary, in most cases at large tan /3, the flipped sign of the amplitude for the decay 
B — » X s ^f is excluded by a combination of the experimental bounds from B s — > n + (JT and 
B s — B s mixing. For low values of tan (3 however, where the constraints provided by 
B s — > /i + fi~ and B s — B s play only a minor role, it is clear that the additional constraint 
provided by B — > X s l + l~ will be fairly useful when attempting to exclude these extreme 
regions of parameter space |33|] . 

Another means of constraining the GFM parameters d^y that has been discussed in the 
literature is via the use of precision electroweak observables, such as mjy and sin 2 9 e g p4|j . 
The constraints in this case typically affect either LL or RR insertions as contributions 
due to LR and RL insertions typically appear at fourth order in the MIA (as opposed to 
second order for LL and RR insertions). Taking into account these corrections, it has been 
shown that, for large values of <5^ L , quite sizeable corrections to both mw and sin 2 e fj 
can be induced that provide a useful limit in constraining extreme values of 5f L [34]. The 
improved measurements that will be made at the LHC and the next linear collider will 
serve to strengthen these constraints even more. 



9. Conclusions 



In conclusion we have discussed the current limits that can be placed on SUSY flavour 
violation. In particular, we have included all the relevant tan /3 enhanced corrections that 
appear at BLO when calculating the SUSY contributions to the various processes under 
consideration. For the decay B — > X s j we have reiterated the need to include BLO 
corrections when deriving limits in the large tan/3 limit, showing the effect of varying 
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parameters such as \i and rrig on the BLO focusing effect discussed previously in 14, 15]. 
For the processes B s — > \i+ \i~ and B s — B s mixing we have shown that useful limits can 
be placed on all four insertions with the current experimental bounds. We have also 
combined the limits from all three processes and investigated the effect the constraints 
have in a variety of scenarios. In particular, we have shown that both B s — > [i+ [x~ and 
B s — B s mixing can play an important role in constraining RL and RR insertions for tan (3 
as low as 30. Our combined analysis also reveals that the case \x < can still be allowed for 
small amounts of SUSY flavour violation, even for light sparticle masses. The constraints 
also serve as a useful tool when constraining multiple sources of flavour violation and 
useful bounds are already available in the 8 ( f jL — 5^ R plane. Finally, we have discussed the 
prospects available for future experiments that aim to probe the processes B s — ► /U + /i~ 
and B s — B s mixing. In particular, we have identified a pattern of possible observations at 
large tan/3 that could only be associated with flavour violation in the RL or RR sectors. 
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